Formation of molecules in an expanding Bose-Einstein condensate 
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A mean field theory of expanding hybrid atom-molecule Bose-Einstein condensates is applied to 
the recent MPI experiments on 87 Rb that demonstrated the formation of ultracold molecules due to 
Feshbach resonance. The subsequent dissociation of the molecules is treated using a non-mean-field 
parametric approximation. The latter method is also used in determining optimal conditions for 
the formation of molecular BEC. 
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Molecular Bose-Einstein condensates (BEC) have been 
recently formed in experiments on atomic BEC 0,0,001 
and on quantum-degenerate Fermi gases These ex- 
periments exploited the effect of Feshbach resonance , 
present when the energy of an atomic pair is close to the 
energy of a bound molecular state. The energy mismatch 
can be controlled by applying an external magnetic field, 
thanks to the difference between the magnetic momenta 
of the molecule and of the atomic pair. By varying the 
magnetic field the energy of the molecular state is forced 
to cross different states of the atomic pairs, belonging 
to a discrete spectrum in the case of a trapped gas, or 
to a continuum otherwise. The resonance value of the 
magnetic field strength Bq corresponds to the crossing of 
the lowest discrete state, or the lower boundary of the 
continuum, respectively. 

In the experiments Q, 0, @ the resonances were 
mostly swept in a backward direction, so that the molec- 
ular state crossed the atomic ones downwards. This led 
to the transfer of population from the lowest atomic state 
in the case of BEC, or from an energy band in the case of 
a Fermigas, to the molecular state, as had been proposed 
in Ref . |7| ■ Assuming all the atomic population is initially 
in the BEC state, the backward sweep would have been 
ideally suitable for forming molecules, were it not for 
two destabilizing mechanisms. The resonant molecule is 
generally populated in an excited rovibrational and elec- 
tronic state, and therefore can be deactivated by inelas- 
tic collisions with atoms and other molecules (see Refs. 
0, El) In addition, during the backward sweep some 
higher-lying non-condensate atomic states can be popu- 
lated temporarily (by counterintuitive transitions see 
Ref. 01) i most notably in a strong resonance or at low 
densities (see Ref. [ll|). These two effects restrict the ef- 
ficiency of conversion from the atomic condensate to the 
molecular one. The non-condensate atoms produced by 
molecular dissociation are formed as entangled atomic 
pairs (see Ref. [u|). This process cannot be described 
by mean field (MF) theories. Therefore the analysis of 
formation of a molecular BEC incorporating both effects 
requires a non-MF approach in which the damping due 
to deactivating collisions is incorporated. Such an ap- 
proach — a generalized parametric approximation (PA) 



— has been described in Ref. and applied to a Na 
BEC. The non-MF effects become even more important 
in forwaed sweeps. An alternative study of molecular 
formation described in Ref. |12|. using the macroscopic 
quantum-dynamics approach of Ref. |13( , takes into ac- 
count non-MF effects but not inelastic collisions. 

In the 87 Rb experiments 0] a BEC of 10 5 atoms has 
been kept initially in a harmonic trap with frequencies 
2-7T x (50, 120, 170) Hz. The magnetic field ramp was 
started, following the switching off of the trap, after a 
pre-expansion interval t p > 2 ms. This measure allows 
the reduction of condensate density, and thus improves 
conversion efficiency (see Ref. H3)- The field has been 
ramped from 1008 G to 1005.2 G with various speeds, 
passing the 1007.4 G resonance of strength A w 0.2 G 
(see Ref. ^ij), and held for 3 ms. During this time a 
gradient magnetic field has been applied, leading to a 
relative motion of the atoms and the molecules due to 
the difference of magnetic momenta between the molecule 
and an atomic pair, // w 2.8/ig (see Ref. |lq). Next, the 
magnetic field was ramped up to 1008 G (in a forward 
sweep), converting the molecules back to atoms. The 
reconverted atoms inherit the velocity of the molecules, 
and thus form a second condensate cloud. 

Association of atoms in an expanding BEC. The PA 
approach as used in [ll| for a homogeneous gas is inap- 
plicable directly to experiments on traps. However, the 
MF approach of (which can be easily extended to 
traps) should suffice for the description of association in 
the case concerned. This is shown by Fig.Q] It compares 
the atomic and molecular densities calculated with the 
PA approach 0] and the MF approach of 0] , for an 
initial atomic density corresponding to the peak density 
reached at the expansion time of 2.3 ms. The values used 
for the deactivation rate coefficients are k a — 7 x 10~ n 
cm 3 /s (see Ref. ^|) an d k m — 5 x 10~ n cm 3 /s for atom- 
molecule and molecule-molecule collisions, respectively, 
The elastic scattering length is a a = 99 atomic units 
(see Ref. 0). As one can see, already when the mag- 
netic field is 0.3 G below the resonance (in a sweep total- 
ing 2.2 G) the temporary occupation of non-condensate 
atom states practically vanishes, and the results of the 
two kinds of calculations coincide. For faster sweeps or 
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molecules are formed after the expansion starts when the 
trap is switched off and the densities decrease substan- 
tially. 

The expansion of a pure atomic condensate cloud has 
been considered in Ref. |l7|. We generalize this theory 
here to the case of a hybrid atom-molecule condensate. 
Let us consider an initial atomic held with a Thomas- 
Fermi distribution, and a zero molecular held, at t < 
to — t p . We can represent the two mean fields in the 
following form, 
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FIG. 1: Fraction of atoms surviving in the atomic conden- 
sate, and those converted to molecules and non-condensate 
atoms, calculated using the PA (dot-dashed and solid lines, re- 
spectively) and the MF approach (dot-dot-dashed and dashed 
lines, respectively) for a homogeneous gas with the initial 
atomic density 3 x 10 13 cm -3 and a magnetic ramp speed of 
1 G/ms. The fraction of non-condensate atoms calculated 
by using the PA is represented by the dotted line. Pluses 
and circles represent the fraction of surviving and converted 
atoms, respectively, calculated using the MF approach but 
taking into account the effects of spatial inhomogeneity and 
expansion. 



higher densities the two results converge even faster. 

The extension of the MF approach of Refs. d, to 
inhomogeneous expanding gases requires the solution of 
two coupled Gross-Pitaevskii equations for the atomic 
ip a (r, t) and the molecular tp m (r, t) mean holds, 



Lp a {r,t) =A(t) $ a (p,t)e 
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(2) 



tp m (r, t) = A (i) $ m (p, t) e 
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using the scaled coordinates pj — rj/bj (t), I < j < 3, 
and a scaling factor A (t) — (b\ (t) &2 (t) &3 (*)) • 
scales bj obey the set of equations 



The 



bj(t)=uj?A 2 (i)/fy(t),&i(*o-tp) 
The phase in Eq. J5J), 
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iUifa (r, t) = 
A-kU 2 



— V 2 + e a (t) + V a (r, t) - -k a \if m (r, t) | 5 



accumulates most of the contributions of the kinetic en- 
ergy. Here eo = AnTi 2 ''a a nQjm is a chemical potential of 
the atomic BEC and no is its peak density while the trap 
is on. 

Substitution of Eq. J2J into leads to the following 
set of coupled equations for the transformed mean fields, 
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-fc m |$ ro (p,i)| ; 



$ m (p,t)+A(t)g<Z> 2 a (p,t). 



Here m is the mass of the atom and e a (t) = 
— |/i (B (t) — Bq) is the time-dependent Zeeman shift of 
the atom relative to half the energy of the molecular 
state. The external magnetic held B (t) is kept constant 
while t < to, and is linearly ramped at t > to so that 
B (0) = i?o- The atoms are considered trapped in a 
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harmonic potential V a 



(T,t) = fE^< 

3=1 



{to — t p — t), 



which is switched off before the ramping, at t = to — t p . 
The atom-molecule coupling is related to the phenomcno- 
logical resonance strength A as \g\ 2 — 2Trh 2 \a a \uA/m 
(see Ref. 0). The molecular trap potential and elastic 
collisions involving molecules can be neglected since the 



The residual kinetic energy terms are neglected here for 
the same reasons as in Ref. ^7|- The loss processes and 
atom-molecule transitions distort <& Q (p, t) and $ m (p, t) 
from the Thomas-Fermi shapes, leading to additional en- 
ergy shifts compared to the pure atomic case. These 
shifts, however, are of the order of A 2 (t) eo, and can only 
lead to a negligible small shift of the resonance (of less 
than 10 _4 G in the present case). 

Equations J5J have a rather clear physical sense. They 
describe a ballistic expansion of the atomic and molecular 
BEC with the same velocity distribution, reducing the 
densities by the factor A 2 (t) , and leading to a rescaling of 
the coupling and deactivation parameters. This reflects 
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FIG. 2: Fraction of atoms surviving in the atomic condensate 
(solid line) and those converted to molecules (dot-dashed line) 
calculated with the MF approach, taking into account spa- 
tial inhomogeneity and expansion, are shown for k a = 10~ 10 
cm 3 /s. The dashed and dotted lines show the fraction of 
surviving and converted atoms, respectively, calculated for 
k a = 7 x 10" 11 cm 3 /s. 



the fact that the acceleration predates the formation of 
molecules that inherit the velocity of the atoms they are 
formed from. 

The results of a numerical solution of Eq. JSJ) are shown 
in Fig. ^ The inhomogeneity and the expansion reduce 
the atomic density, leading to a slower loss of atoms and 
molecules. The dependence on the sweep rate is pre- 
sented in Fig. [21 for two values of k a (the best fit and 
upper limit taken from Ref. 0). The results are m 
agreement with the experimental data of Ref. Q report- 
ing that ~ 7% of atoms are converted to molecules and 
~ 30% remain in the atomic BEC for ramp speeds less 
than 2 G/ms. 

Dissociation of molecules on a forward sweep. Con- 
sider now the dissociation of molecules on a forward 
sweep, when the molecular state crosses the atomic ones 
upwards. Such a sweep has been used in experiments 
for detection of molecules. The experiments demon- 
strate no significant dependence of the reconversion ef- 
ficiency on the ramp speed. Our calculations using the 
PA [llj demonstrate that under the experimental con- 
ditions all molecules are dissociated except for a small 
part lost due to deactivating collisions during the sweep. 
The molecules dissociate to en tang led atomic pairs in a 
wide energy spectrum (see Ref. [llj) characterized by the 
peak energy E pea k. The calculations demonstrate that 
Epeak increases with the ramp speed (see Fig. yj. The 
results show no significant dependence on the molecular 
density, justifying the applicability of the homogeneous 
density approximation to real inhomogeneous situations. 
Although the ramp-speed dependence of -Epoak is clearly 
nonlinear, within the interval 0.2 G/ms < B < 5 G/ms 
it can be approximated by a straight line with the slope 
of 70 nK ms/G, which is close to the experimental value 
of 60 nK ms/G reported in Ref. (lif . 

In the present case the dissociation can be described by 



FIG. 3: Peak energy of the non-condensate atoms formed by 
the dissociation of molecules in a forward sweep calculated 
using Eq. @ (dashed line) and the PA, for the initial molec- 
ular densities 10 10 (o), 10 11 (x) and 10 12 (+) cm" 3 . The solid 
line represents the results of a linear fit. 



a simple analytical curve-crossing model based on the ap- 
proach of Ref. . A molecule dissociates into two entan- 
gled atoms, each with an energy E, by a crossing occur- 
ring at /i (B (t) — Bo) = 2E. In the weak-resonance case, 
the quantum curve-crossing theory and the semi- 
classical Landau-Zener theory yield approximately the 
same result for the crossing probability, i.e., e 27rA -1« 
l_ e -2*\ ^ 27r ^ w here A = %irfia a n m (t) A/ (m\B\j (see 

Ref. Neglecting deactivating collisions the molecu- 

lar density n m (t) = \ip m (t) | 2 is given by the solution of 
Eq. (50) in Ref. 



n m (*) = n m (td) exp 



t d 



where td is the starting time of the dissociation ramp. 
The resulting energy distribution of the produced atoms 
is given by 



f(E) = lVEE pc 3 i 2 exp 



~3 ( E / E v Qak ^ 3/2 



with 



E, 



peak 



hm\B\ 
4a„A 



2/3 



(6) 



(7) 



This model is in good agreement with the numerical re- 
sults (see Fig. |3J). Similar models have been indepen- 
dently developed in Refs. P-H^]. 

Optimized molecular formation on a backward sweep. 
Figure ^ demonstrates that the molecular population 
reaches a maximum, corresponding to a conversion ef- 
ficiency ~ 50%, about 0.1 G below the resonance, and 
falls by half about 0.5 G below the resonance. The low 
conversion efficiency observed in experiments Q is due 
to collision loss during the long sweep to 2.2 G below the 
resonance. Although the PA results of Fig. ^ overesti- 
mate the loss by neglecting the expansion, the molecular 
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FIG. 4: Conversion efficiency (a) and the lifetime of the molec- 
ular BEC T m (b) at the optimal ramp speed (dB/dt) t , all 
plotted as functions of the initial atomic density no, calculated 
for the resonances at 1007 G (circles) and 685 G (triangles), 
using the rates of molecule- molecule deactivation 10 -10 (solid 
lines) and 0.5 x 10 -10 cm 3 /s (dashed lines). The dash-dotted 
and dotted lines in part (b) demonstrate (dB/dt) t for the 
two rates, respectively. 




FIG. 5: Conversion efficiency as a function of the ramp speed 
dB/dt, calculated at the initial atomic density 10 13 cm -3 for 
the resonances at 1007 G (circles) and 685 G (triangles), using 
the rates of molecule-molecule deactivation 10 -10 (solid lines) 
and 0.5 x 10~ 10 cm 3 /s (dashed lines). 



density and the resonance strength (see Fig. 0J. 

Conclusions. The experiments Q on formation of 
molecules from a 87 Rb BEC due to a Feshbach reso- 
nance can be described by a MF theory of expanding 
atom-molecule BEC. Non-MF calculations using the PA 
describe the energy of non-condensate atoms formed by 
molecular dissociation in a forward sweep and can pre- 
dict optimal conditions for the formation of a molecular 
BEC in a backward sweep. 



population of the expanding gas falls by half on reaching 
2.2 G. Under the conditions of the experiments the life- 
time of the molecules formed is about 0.5 ms, too short 
for an effective detection. As in the case of Na (see llj), 
the lifetime and conversion efficiency increase on reduc- 
ing the initial condensate density. This is demonstrated 
in Fig. 01 showing results of calculations using the PA 
fll| for the homogeneous case. The use of this non-MF 
approach is important for such situations, as the MF ap- 
proach becomes inadequate in the vicinity of the peak 
molecular occupation (see Fig.^l, and at low condensate 
densities. A use of the weaker resonance at 685 G with 
A w 17 mG /it = 1.4/is (see Ref. 15]) should increase the 
conversion efficiency (see Fig.0J. 

As in the case of Na studied in Ref. , the conver- 
sion efficiency is determined by a concurrence of three 
processes: the association of the atomic BEC and the 
two loss processes — the dissociation of the molecular 
BEC onto non-condensate atoms and the deactivation 
by inelastic collisions. A reduction of the ramp speed 
enhances all three processes. The calculations of Ref. 
[l2|. taking no account of the inelastic collisions, result 
in a monotonic increase of the conversion efficiency with 
a decrease of the ramp speed. The inelastic collisions, 
however, tend to reduce the conversion efficiency at slow 
ramp speeds (see Fig.|3J), as has also been demonstrated 
in Ref. [ll| . The optimal ramp speed increases with the 
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